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ABSTRACT ^_^ 

In recent work ( [Hayward et al.|20lTa) we have suggested that the high-redshift (z ~ 2 — 4) 
bright submillimetre galaxy (SMG) population is heterogeneous, with major mergers con- 
tributing both at early stages, where quiescently star-forming discs are blended into one 
submm source ('galaxy-pair SMGs'), and late stages, where mutual tidal torques drive gas 
inflows and cause strong starbursts. Here we combine hydrodynamic simulations of major 
mergers with 3-D dust radiative transfer calculations to determine observational diagnos- 
tics that can distinguish between quiescently star-forming SMGs and starburst SMGs via 
integrated data alone. We fit the far-IR SEDs of the simulated galaxies with the optically 
thin single-temperature modified blackbody, the full form of the single-temperature modified 
blackbody, and a power-law temperature-distribution model. The effective dust temperature, 
Td, and power-law index of the dust emissivity in the far-IR, /?, derived can significantly de- 
pend on the fitting form used, and the intrinsic f3 of the dust is not recovered. However, for 
all forms used here, there is Xd above which almost all simulated galaxies are starbursts, so a 
Td cut is very effective at selecting starbursts. Simulated merger-induced starbursts also have 
higher LiR/M gas and Ltr/Lfuv man quiescently star- forming galaxies and lie above the 
star formation rate-stellar mass relation. These diagnostics can be used to test our claim that 
the SMG population is heterogeneous and to observationally determine what star formation 
mode dominates a given galaxy population. We comment on applicability of these diagnostics 
to ULIRGs that would not be selected as SMGs. These 'hot-dust ULIRGs' are typically star- 
burst galaxies lower in mass than SMGs, but they can also simply be SMGs observed from a 
different viewing angle. 

Key words: dust, extinction - galaxies: high-redshift - galaxies: starburst - infrared: galaxies 
- radiative transfer - stars: formation. 



1 INTRODUCTION 

1.1 The two modes of star formation 

Star formation is one of the fundamental processes driving galaxy 
formation: it depletes the gas content of galaxies, enriches the 
interstellar medium (ISM) with metals, and deposits energy and 
momentum via supernovae, stellar winds, and radiation pressure. 
Furthermore, the light emitted by stars encodes much information 
about the current physical properties of a galaxy and the galaxy's 
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formation history. Thus understanding star formation is crucial for 
understanding galaxy formation and evolution. 



An important step toward understanding the star formation 
processes that built up galaxies across cosmic time is determin- 
ing where and when most stars are formed, be it in disc galaxies 
or in starbursts triggered by, e.g., galaxy mergers, which are short- 
lived but can dramatically alter a galaxy's properties. An increasing 
amount of observational evidence supports the notion that there are 
two modes of star formation (e.g., Wuyts et al. 201 lab, Elbaz et al. 
|20l71|Rodighiero et al.|201 l[|Nordon et al.|2012| >, typically referred 
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to as quiescently star-forming or quiescen{j(that occurring in nor- 
mal disc galaxies) and starburst (found in unstable discs and merg- 
ing galaxies at first passage and coalescence, though whether a star- 
burst is induced in the latter depends on factors such as gas content, 
orbit, and mass ratio of the progenitors; e.g., |Cox et al.|2 008). One 
clear difference is that gas depletion time- scales of starbursts are 
significantly lower than those of quiescently star-forming galaxies. 
Star formation in starbursts tends to be dominated by the nuclear 
regions, whereas quiescent star formation is more extended. Addi- 
tionally, starbursts may obey a different global Kennicutt-Schmidt 
(KS) relation ( [Kenni cutt 1998a, Schmidtp959j than quiescently 
star-forming disc galaxies: Daddi et al. (2010) and Genzel et al. 
(2010) argue that the normalisation of the KS relation for starbursts 
is ~ 4 — 10 x greater than that for quiescently star-forming discs. 
However, this conclusion strongly depends on the bimodality of the 
values adopted for the CO-H2 conversion factor, so the data in fact 
may be consistent with a single KS law normalisation ( Narayanan 
|etal.|20lTc] >. The ratio of infrared (IR) luminosity to molecular gas 
mass, which is a measure of the global star formation efficiency of 
the system, is larger in starbursts than in quiescently star-forming 
galaxies, though the magnitude of the difference is also sensitive 
to the CO-H2 conversion factor. Furthermore, the relationship be- 
tween SFR and dust mass may also show a bimodal behaviour (da] 
|Cunha et al.|2010) . 

At a given redshift, most galaxies lie on a tight relation be- 
tween SFR and stellar mass (M*; |Noeske et aT]|2007a|b| |Daddi| 
|et al.|2007||Rodighiero et al.|2010||Karim et al.|2011| >. The relation 
arises because star formation is supply-limited, so, on average, SFR 
correlates well with cosmological gas accretion rates, which are 
well-correlated with halo mass ( Keres et al. 2005 2009 ; Faucher- 
|Giguere et al.|20li] ). In this picture, starbursts are transient events 
that cause a galaxy to move significantly above the SFR-M* rela- 
tion for a short (~ 50 — 100 Myr) time. During the burst the gas is 
rapidly consumed, the SFR declines sharply, and the galaxy returns 
to the SFR-M* relation or is quenched (i.e., drops significantly be- 
low the relation), depending on factors such as merger mass ratio, 
orbit, and gas fraction. 

Though there is some observational support for two star for- 
mation modes, the underlying physics is not fully understood. Thus 
further detailed observations are crucial. However, it can be dif- 
ficult to observationally determine which mode of star formation 
dominates a given galaxy population; this complicates efforts to 
understand the underlying physics. This is especially a problem 
at high redshift because of stricter observational limitations, and 
lessons learned from the local universe may not apply to high- 
redshift galaxies. For example, at high redshift gas accretion rates 
are significantly higher than locally (e.g.JKeres et al.|2005|), so gas 
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fractions (Erb et al. 2006, Tacconi et al. 2006 2010, Daddi et al. 
2010) and star formation rates (Noeske et al. 2007a b, Daddi et al. 



2007 ) of galaxies at fixed galaxy mass increase rapidly with red- 
shift. Consequently, even a typical star-forming galaxy at z ~ 2 
can reach ULIRG lumin osities (e.g., [Dadd i et al. 2005 2007l[Hop^l 
|kins et al. 12008] [20T0l |Dannerbauer et al.|2009) . It would thus be 
useful to have simple observational diagnostics that can be used to 
determine which mode of star formation dominates a given galaxy 
or galaxy population. This is one of the goals of this paper. 



1 Confusingly, the term 'quiescent' is also used to refer to galaxies that 
have essentially no ongoing star formation; here the term 'quiescent' always 
means 'quiescently star-forming'. 



Submillimetre galaxies (SMGs; 



Smail, Ivison, & Blain| [T997] 



Barg er et al.|[T998| |Hughes et al.||1998| |Eales et al.||1999T "see 



Blain et al. 2002 for a review) are a class of high-redshift (me- 
dian z ~ 2.3 - 2.6; |Chapman et~aL] [20051 |Yun et al.||25T2] > 
galaxies notable for their extreme lum inosities (bolometric lumi- 
nosity L bo i ~ 10 12 - 10 13 L ; e.g., |Kovacs et al.||2006| |Mag-| 
|nelli et al.|2010| |2012| ), almost all of which is emitted in the IR. 
Since they seem to be powered by star formation rather than ac 



tive galactic nuclei (AGN; [Alexander et al.|200 5a b , 2008 , Valiante 



et al.||2007| [M enendez-Delm estre et al.||2007| |2009[ |Pope et~ 



2008) [Younger et al.||2008| |2009c), they have inf erred SFRs of 
10 2 - 10 a M yr x (e.g., |Coppin et al.|2008| ), much greater 



than those of even the most extreme local galaxies. 

It has long been known that during major mergers tidal torques 
drive significant amounts of gas inward, resulting in a strong burst 
of star formation pernquist|1989[|Barnes & Hernquist|1991|[T996l 
Mih os & Hernquist||1996| >. Locally, ultra-luminous IR galaxies 
(ULIRGs, defined by Lir > 10 12 L ) are exclusively late-stage 
major mergers (e.g., Sanders & Mirabel 1996, Veilleux, Kim, & 
Sanders 2002, Lonsdale, Farrah, & Smith 2006), so the strong gas 
inflows induced during the coalescence stage of a major merger 
seem necessary to power the most luminous and rapidly star- 
forming galaxies. 

The identification of ULIRGs with late-stage major mergers 
in the local universe has caused many researchers to suspect that 
SMGs, some of the most IR-luminous galaxies at high redshift, 
are also late-stage major mergers. There is much observational evi- 
dence that supports this picture (e.g., Ivison et al. 2002 , 2007 2010 , 
Chapman et al. 2003]|Neri et al.|20031|Smail et al.|2004[|Swinbank 



etal. 2004, Greve et al.|2005 Tacconi et all2006.. 2008, Bouche 



eTaTl|2007l |Biggs & Ivison| [20081 [Capak et al.|[2Q08| [Younger 
eT"aLp0()8l |2010| |Iono et al. ||2009[ |Engel et al.poTof |Riechers 
et al.|2011a[btTBussmann et al.||2012||Magnelii et al.|2012). Fur- 



thermore, by combining hydrodynamic simulations with dust ra- 
diative transfer (RT), we have shown that simulated major merg- 
ers have observed 850-|jjtl fluxes and typical spectral energy dis- 
tributi ons (SEDs;|Narayanan et al.||2010b[ |Hayward et al.|2011a[ 
but cf. Chakrabarti et al. 2008 ), stellar masses (Michalowski et al. 
2011 ), and CO properties (Narayanan et al. 2009) consistent with 
observed SMGs. Semi-analytic models (SAMs) typically also find 
that merger-induced starbursts (though not necessarily major merg- 
ers, as minor-merger-induced starbursts dominate in some models) 
account for the bulk of the SMG population {[Baugh et al.|[2005 



Fontanot et al. 2007, Swinbank et al. 2008, Lo Faro et al. 2009 
Fontanot & Monaco 2010, Gonzalez et al. 2011, but cf. Granato 
et al. 2004 ). However, Dave et al. ( 2010 ) have claimed that there are 
not enough major mergers to account for the observed SMG popu- 
lation. Instead, they argue that a significant fraction of the popula- 
tion must be massive discs fuelled by smooth accretion and minor 
mergers. 

In Hayward et al. (201 la hereafter HI 1) we suggested a mod- 
ification to the canonical picture, arguing that SMGs are not purely 
late- stage major mergers but rather a heterogeneous population, 
composed of: 1. late-stage major mergers undergoing strong star- 
bursts; 2. early- stage major mergers ('galaxy-pair SMGs') which 
are quiescently forming stars; 3. physically unrelated galaxies 
blended into one submm source ( [Wang et"aT]|2011) ); and 4. iso- 
lated disc galaxies and minor mergers, which only contribute sig- 
nificantly at the fainter end of the population. One reason early- 
stage mergers also contribute is that observed submm flux increases 
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Figure 1. The first three rows show mock observed-frame 850-fim continuum images (assuming z = 2) of the merger from Hll near apocentre (first row; 
time t = 0.28 Gyr; nuclear separation c^bh = 71/i _1 kpc; quiescently star- forming), during final approach (second row; t = 0.63 Gyr; c^bh = 13/i _1 
kpc; quiescently star-forming), and at the peak of the starburst (third row; t = 0.7 Gyr; c?bh = 0.1/i -1 kpc; starburst). The fourth row shows the isolated 
disc at t = 0.28 Gyr (quiescently star-forming). The simulated images have been convolved with a Gaussian filter in order to mimic the intrinsic simulation 
resolution (first column; FWHM = 0.2/i -1 kpc) and the resolution achievable with various telescopes/interferometers: The second column has FWHM = 0.5 
arcsec (3/i -1 kpc at z = 2), the resolution roughly characteristic of (sub)mm interferometers such as the SMA, the IRAM Plateau de Bure Interferometer, 
and ALMA. (Note we have not made any attempt to model the interferometric process.) See Fig.[2]for zoom-ins of the boxed regions. Third column: FWHM = 
7 arcsec (42/i _1 kpc), representative of the resolution available with, e.g., Herschel PACS and the current configuration of the LMT. Fourth column: FWHM 
= 15 arcsec (89/i -1 kpc), representative of, e.g., Herschel SPIRE and the JCMT. The fields of view are 23.5 arcsec (140/i -1 kpc) for the first three columns 
and 59 arcsec (351/i -1 kpc) for the last. For all but the first column, Gaussian noise with realistic amplitude has been added. The same noise map is used for 
all images in a given column. When present, scale bars show the beam FWHM, and contours correspond to 3a, 4<r, 5cr, 10cr, and 20cr. The first column has 
logarithmic scaling; all others are linear. Units are arbitrary. 



© 2012 RAS, MNRAS 000,[l]J2T] 



4 Hayward et al 




Figure 2. Zoom-ins of the regions marked with white boxes in the second 
column of Fig.[T] The upper left (right) panel shows the left (right) region 
from the first-row image of Fig. [T] The middle left corresponds to the second 
row, the middle right to the third row, and the bottom to the fourth row of 
n The FOV of each image is 3 arcsec (18/i _1 kpc). See the caption of 
If or further details. 



row, int 

Fig.[T] r 

Fig.p 



rather weakly with SFR, and the starburst mode is significantly less 
efficient at boosting submm flux than the quiescent mode. Phys- 
ically, submm flux scales more weakly in starbursts for two rea- 
sons: 1. In high-redshift mergers significant star formation occurs 
before the starburst induced at the coalescence of the two discs. 
During the starburst the 'old' stars formed pre-burst can provide 
a non-negligible contribution to Lir. For example, for the merger 
shown in the right panel of fig. 1 of HI 1, the stars formed pre-burst 
account for ~ 1/7 of the total IR luminosity at the time of the 
burst. Consequently, Lir does not scale linearly with the instan- 
taneous SFR during the burst because of the dust heating from the 
older stars (see, e.g., section 2.5 of Kennicutt 1998b and section 3.1 
of Hll). 2. Driven primarily by the strong drop in dust mass dur- 
ing the merger, the effective dust temperature of the SED increases 
sharply during the starburst. This effect mitigates the increase in 
submm flux caused by the increased luminosity. These two effects 
result in a very weak scaling (Ssso oc SFR 01 , compared to Ssso oc 
SFR 4 for the quiescently star-forming galaxies). For the example 
in Hll, an increase in SFR of ~ 16 x causes a < 2x boost in the 
submm flux. 

The other reason early- stage mergers contribute significantly 
is that the large (FWHM ~ 15 arcsec, or ~ 120 - 130 kpc at 



z rsj 2 — 3) beams of the single-dish submm telescopes used to 
detect SMGs cause the two merging discs to be blended into a sin- 
gle source for much of the pre-coalescence stage of the merger. 
Fig. [T] demonstrates the effects of blending by showing observed- 
frame 850- |jjtl (assuming z = 2) continuum images of the merger 
from Hll near apocentre (first row), during final approach (second 
row), and at the peak of the starburst (third row). The fourth row 
shows the isolated disc from Hll. The first column shows the full- 
resolution simulated images. The others are meant to mimic the 
resolution attainable with various telescopes/interferometers. See 
the figure caption for full details. 

Comparison of the columns shows the importance of blend- 
ing: For instruments with resolution ~ 15 arcsec (fourth col- 
umn of Figure [TJ, at all stages (even at apocentre) the two discs 
are blended into one submm source, and the morphologies of the 
galaxy-pair SMGs, the late- stage merger-induced starbursts, and 
the isolated disc are similar. This resolution is typical of the sur- 
veys that have been done to date with SCUBA, AzTEC, and Her- 
schel SPIRE, so those surveys should contain many galaxy-pair 
SMGs consisting of widely separated merging discs blended into 
one submm source. With resolution ~ 7 arcsec (third column of 
Figure [TJ the widely separated galaxy-pair SMGs (first row) can be 
resolved, but some mergers which are nearer coalescence but still 
predominantly forming stars quiescently will be blended into one 
source (second row). Thus wide-field surveys done with, e.g., Her- 
schel PACS or the Large Millimeter Telescope (LMT) can more 
effectively distinguish galaxy-pair SMGs from starburst SMGs and 
should have a higher fraction of close pairs in the FIR/(sub)mm 
than the SCUBA, AzTEC, and Herschel SPIRE surveys. When the 
resolution is ~ 0.5 arcsec (second column of Figure [T] see Fig. 
p] for zoom-ins of the boxed regions), typical of interferometers 
such as the SMA, the IRAM Plateau de Bure Interferometer, and 
ALMA, all but the mergers nearest coalescence are resolved into 
two components, as this beam size corresponds to ~ 3/i _1 kpc 
at z ~ 2 — 3. At such low separations the mergers typically are 
undergoing strong starbursts, so the SMGs observed as single com- 
ponents with interferometers should be predominantly starbursts or 
isolated discs. Note that when widely separated merging discs are 
observed with some interferometers one of the sources may be out- 
side the FOV, so from the interferometry such SMGs would appear 
to be isolated discs whereas they are actually galaxy-pair SMGs. 
Clearly it is difficult to disentangle the SMG sub-populations using 
the number of components observed or the morphology. 

The blending of the discs during the early stages of a merger is 
effective at creating SMGs because, by treating the system as a sin- 
gle source, the integrated submm flux and SFR are both doubled; 
see sections 3.2 and 4.1 of Hll for further details. This blending, 
along with the stronger scaling between submm flux and SFR in the 
quiescent mode compared to the starburst mode, causes early- stage 
mergers to contribute significantly to the SMG population (details 
of the contribution will be presented in Hayward et al., in prepa- 
ration). The brightest SMGs should still be merger-induced star- 
bursts, but early-stage mergers must provide a significant contribu- 
tion to the population. 

There is already much evidence that some SMGs are early- 
stage mergers. For example, [Engel et aT1j2010) used submm inter- 
ferometry to show that approximately half of their SMG sample are 
well-resolved binary systems. (See also Tacconi et al.|2006 2008 , 



|Bothwell etaL 2010, Rie chers et al.|2011a|b| ) Two of the twelve 
SMGs in the Engel et al. sample consist of two well- separated, re- 
solved components (projected separations > 20 kpc), and it is pos- 
sible that they have missed one component of galaxy-pair SMGs 
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with more widely separated components because of the limited 
field of view. In such widely separated systems the star formation 
induced by the tidal torques exerted by the discs upon one another 
is not sufficient to drive a strong starburst, so the discs would form 
stars at similar rates even if the companion were absent. Thus such 
systems should be considered physically analogous to normal disc 
galaxies rather than late- stage mergers because they are forming 
stars via the quiescent rather than starburst mode. In addition to 
the evidence from CO interferometry, support for the galaxy-pair 
contribution is provided by the frequency of multiple counterparts 
in the radio (e. g.,|Ivison e t al. 2002 2007 ; Chapman et al. 2005 , 
|Clements et al.|2008| | Younger 'et alj2009a||Yun et al.|2012|> and 
24-iim (e.g., |Pope et al.|20061|Hwang et al.|2010[|Yun et al.|2012) 
emission and SMGs with morphologies that appear more like discs 
than late-stage mergers (e.g., Bothwell et al. 2010] |Carilli et al.| 
[20T0l[Ricc7ardelli et al.|2010[|Targett et al.|201l| T(Note, however, 
that in gas-rich mergers discs can rapidly re-form, potentially con- 
fusing interpretation of these results; Springel & Hernquist 2005, 
Nara yanan et al.|2008b[ [Robertson & Bullock 2008[ |Hopkins et al.| 



20091 ) 



We have shown that physical arguments and observations sug- 
gest that the SMG population is a mix of quiescently star-forming 
galaxies and starbursts. Thus SMGs differ significantly from local 
ULIRGs, which are exclusively starburst- or AGN-dominated, so 
one should draw comparisons between the two populations with 
care. The heterogeneity complicates physical interpretation of the 
SMG population. For example, one should not apply a CO-H 2 
conversion factor appropriate for starbursts to the quiescently star- 
forming sub-population of SMGs. Furthermore, proper treatment 
of all sub-populations is key for reproducing the observed SMG 
number counts (Hayward et al. 2011b Hayward et al., in prepa- 
ration). However, the relative contributions of these various sub- 
populations is not yet observationally well-determined, and predic- 
tions of the relative contributions depend sensitively on uncertain 
model details. We do not predict the relative contributions in the 
present work. Instead, we wish to determine how one can obser- 
vationally distinguish between starburst-driven (late- stage merger) 
SMGs and those powered by quiescent star formation even when 
only integrated data are available. One can then use these diagnos- 
tics to observationally constrain the relative contributions of star- 
bursts and quiescently star-forming galaxies to the SMG popula- 
tion, thereby testing the bimodality we claim exists. 

The rest of the paper is organised as follows: We describe 
our simulation methodology in Section [2] Section [3] presents mul- 
tiple observational diagnostics that can be used to distinguish be- 
tween quiescently star-forming galaxies and starbursts from inte- 
grated data alone, including the luminosity-effective dust tempera- 
ture relation (Sect ion | 3.1| ), star formation efficiency (Secti on 1 3. 2) , 
IR excess (Section [^JTand the SFR-M* relation (Section |34]Tln 
SectionHlwe discuss some implications of our work, and in Section 
|5]we summarise and conclude. 



(SPH) simulations of mergers of equal-mass disc galaxies with the 
SUNRISE ( |Jonsson|20"06||Jonsson, Groves, & Cox^OlO] ) polychro- 
matic Monte Carlo dust RT code to calculate synthetic SEDs of 
the simulated galaxies. The combination of Gadget- 2 and SUN- 
RISE has been used to successfully reproduce the SEDs/colours of a 
variety of galaxy populations, both low- and high-redshift, includ- 
ing: local SINGS (Kenni cutTet al.|2003[|Dale et al.|2007|) galax - 

ies (Jonsson et al. 120101); local ULIRGs ([Younger et al.||2009b}; 

» » ■ t „ . » » * 

extended UV discs (Bush et al. 2010); 24 ^mi-selected galaxies 

(Narayana n et al.||2010a|); massive, quiescent, compact z ~ 2 
galaxies (Wuyts et al. 2009 2010); and K+ A/post- starburst galax- 
ies (Snyder et al. 201 1 ), among other populations. These successes 
support our application of GADGET and SUNRISE to modelling 
high-redshift ULIRGs. 



2.1 Hydrodynamic simulations 

Gadget- 20 (Sprin gel et aT]|200l[ | Springel |[2005] > is a TreeSPH 
(Hernquist & Katz 1989) code that computes gravitational inter- 
actions via a hierarchical tree method (Barnes & Hut 1986) and 
gas dynamics via SPH (Lucy 1977, Gingold & Monaghan 1977, 
Springel 2010). The formulation of SPH used explicitly conserves 
both energy and entropy when appropriate (Sprin geT& Hernquist| 
2002 ). Radiative heating and cooling is included in GADGET-2 fol- 



lowing Katz, Weinberg, & Hernquist ( 1996 ). Star formation is im- 
plemented using a volume-density-dependent KS law (Kennicutt 
Il998al, psfr oc pg a ^, with a minimum density threshold. The as- 
sumed KS index results in a surface-density-dependent KS law con- 
sistent with observations of z ~ 2 discs ( Krumholz & Thompson 
2007 , Narayanan et al. 2008a 201 la), suggesting that it is reason- 
able to use this prescription in our simulations of z ~ 2 mergers. 
The gas is enriched with metals assuming each particle behaves 
as a closed box, so those gas particles with higher SFRs are more 
rapidly metal-enriched. 

We use the sub-resolution two-phase ISM model of Springel 
& Hernquist (2003). In this model, cold, dense clouds are em- 
bedded in a diffuse, hot medium. Supernova feedback (Cox et aL] 
2006 ), radiative heating and cooling, and star formation control the 
exchange of energy and mass in the two phases. A simple model for 
black hole (BH) accretion and AGN feedback (Springel, Di Mat- 
|teo, & Hernquist||2005l |Di Matteo, Springel, & Hernquist||2005] > 
is included. BH sink particles with initial mass 10 5 M are in- 
cluded in both initial disc galaxies. They accrete via Eddington- 
limited Bondi-Hoyle accretion (Hoyle & Lyttleton 1939, Bondi 
& Hoyle 1944). The luminosity of each BH is calculated from 
the accretion rate Mbh assuming the radiative efficiency appro- 
priate for a Shakura & Sunyaev (1973) thin disc, 10 per cent. Thus 
£boi = O.lMc 2 . Five per cent of the luminosity emitted by the 
BHs is deposited into the surrounding ISM. 



2 SIMULATION METHODOLOGY 



We analyse G adget-2 ( |Springel, Yoshida, & White| [2001 , 
Springel 2005) 3-D N-body/smoothed-particle hydrodynamic^] 



2 Recent wo rks JAgertz et a l. 2007 , Springel 2010 Bauer & Springe i|20lT| 
Sijacki et al. 2011) suggest that the SPH method has several significant 
flaws, including artificially suppressing fluid instabilities, preventing effi- 
cient gas stripping of inf ailing structures, and artificially damping turbulent 



eddies in the subsonic regime. Comparison of cosmological simulations run 
with Gadget- 2 to otherwise identical simulations run with the more ac- 
curate moving-mesh code Arepo (Springel 2010) shows that significant 
differences exist even though all physics incorporated in the simulations is 
identical ( Vogelsberge r et al.|20lT)|Keres et al.|2011[[Torrey et al.|2011) . 
However, preliminary comparison of idealised merger simulations of the 
sort presented here run with Gad get- 2 and Arepo suggests that the two 
methods yield similar results (Hayward et al., in preparation). Thus we ex- 
pect our results to be robust in this regard. 

3 A public version of Gadget-2 is available at http: //www. 
mpa-garching .mpg . de/ ~volker/gadget /index . html 



© 2012 RAS, MNRAS 000,[l]J2T] 



6 Hayward et al 



The simulations are initialised in the following manner: Expo- 
nential discs with initial gas fraction f g = M gas /(M* + M g as ) = 
O.qj are embedded in dark matter haloes described by a |Hern- 
|quist| ( |1990| > profile. The progenitor discs are scaled to z ~ 3 as 
described in Robertson et al. (2006) so that the mergers occur at 
z ~ 2. We have selected galaxy masses representative of the SMG 
population (e.g., Michalowski et al. 201 1 ). The gravitational soften- 
ing lengths are 200 h~ pc for the dark matter particles and \00h~ 1 
pc for the stellar, gas, and BH particles. Each disc galaxy is com- 
posed of 6 x 10 4 dark matter particles, 4 x 10 4 stellar particles, 
4 x 10 4 gas particles, and 1 BH particle. Two identical discs are ini- 
tialised on parabolic orbits with initial separation Ri n a = 5R v [ r /S 
and pericentric distance twice the disc scale length (Robertson et al. 
2006). We analyse only equal-mass mergers because these simula- 
tions provide sufficient examples of quiescently star-forming galax- 
ies (during the early stages) and starbursts (near coalescence). The 
differences between quiescent and starburst modes are insensitive 
to orbit and merger mass ratio. (Not all mergers induce strong star- 
bursts, but those are irrelevant for our present purposes.) The phys- 
ical parameters of the simulated major mergers are summarised in 
Tabled] 



2.2 Radiative transfer 



We use the 3-D Monte Carlo dust RT code SUNR!S^1(Jo nsson| 
2006, | Jonsson et al.|201Q) in post-processing to calculate the far- 
UV-mm SEDs of each simulated merger at 10 Myr intervals. We 
briefly describe the Sunrise calculation here, but the reader is en- 
couraged to see Jonsson et al. (2010) for full details. SUNRISE uses 
the stellar and BH particles from the GADGET-2 simulations as ra- 
diation sources. Each stellar particle with age > 10 Myr is treated 
as a single-age stellar population and assigned a STARBURST99 
(Leith erer et al.|[T999) SED template appropriate for its age and 
metallicity. Stellar particles with age < 10 Myr are assigned a tem- 
plate from Gro ves et~aL] {2008), which includes emission from the 
HII regions surrounding the clusters. We do not include the photo- 
dissociation regions for the reasons discussed in detail in section 
2.2.1 of Hll. The stars in the initial discs are assigned ages by as- 
suming that the population was formed at a constant rate equal to 
the SFR of the initial snapshot. The gas and stars present in the 
initial discs have metallicity Z = 0.015, which results in the galax- 
ies being roughly on the z ~ 2 mass-metallicity relation during the 
starburst. The BH particles are assigned SEDs using the luminosity- 
dependent templates of Hopki ns et al.|p007| ). These templates are 
derived from observations of un-reddened quasars, so they include 
the intrinsic power-law emission and reprocessed hot dust emission 
from the torus. 

SUNRISE calculates the dust distribution by projecting the 
Gad GET- 2 gas-phase metal density onto a 3-D adaptive-mesh- 
refinement grid using the SPH smoothing kernel and assuming a 
dust-to-metal density ratio of 0.4 ( |Dwek| 1998) | James et al.|2002) . 
A minimum cell size of 55/i -1 pc is used; this is sufficient to ensure 
the SEDs are converged to within < 10 per cent at all wavelengths. 
Grain compositions, size distributions, and optical properties are 
given by the Milky Way R = 3.1 dust model of Weingartn er &| 



Draine (2001 ) as updated by Draine & Li (2007). The opacity curve 
for this model has a power-law slope in the far-IR f3 w 2. 

To perform the RT we use 10 7 photon packets for each stage, 
~ 10 x the number of grid cells. This limits Monte Carlo noise to 
less than a few percent. SUNRISE randomly emits the photon pack- 
ets from the sources and randomly draws interaction optical depths 
using the appropriate probability distributions. At the interaction 
optical depth a fraction of the photon packet's intensity is absorbed; 
the remainder is scattered into a direction randomly drawn using the 
scattering phase function. This is repeated until the photon packet 
leaves the grid or its intensity drops below a minimum value. 

The energy absorbed by the dust is re-radiated in the IR. SUN- 
RISE assumes a ll dust (except for h alf of the PAHs with grain size 
< 100 A; see Jonsson et al.|2010 for details) is in thermal equi- 



described below. 



Note, however, that we discard all snapshots with f g < 0.4 for reasons 



G SUNRISE is publicly available at http : / /code . google . com/p/ 
sunrise/ 



librium, so the dust temperature is calculated by setting the lumi- 
nosity absorbed by each grain equal to the energy emitted by the 
grain. The equilibrium temperature of a grain depends on the local 
radiation field heating the grain and its absorption cross section, so 
there are in principle n ce iis * %rain sizes different dust temperatures 
in a given SUNRISE calculation. This is important to keep in mind 
when one considers fitting IR SEDs with modified blackbodies, as 
discussed below. 

In high-density environments the ISM can be optically thick 
in the IR; this is common in the central regions of the late-stage 
mergers modelled here. Consequently, one must account for atten- 
uation of the dust re-emission (aka dust self- absorption). Further- 
more, since the IR emission absorbed heats the dust, one must it- 
erate the dust temperature calculation and RT of the dust emission 
until the dust temperatures are converged. SUNRISE uses a refer- 
ence field technique similar to that of Juvela (2005) to perform this 
iteration. We encourage the interested reader to see | Jons son et al.| 
J2QT0] > and |Jonsson & Primack| ( [20T0l ) for details. 

The SUNRISE calculation yields spatially resolved, multi- 
wavelength (for these simulations there are 120 wavelengths sam- 
pling the UV-mm range) SEDs for each galaxy snapshot observed 
from 7 different viewing angles distributed uniformly in solid an- 
gle. The data are analogous to data yielded by integrated field unit 
(IFU) spectrographs. For this work we spatially integrate to calcu- 
late integrated SEDs for the system. When calculating observed 
flux densities we assume the simulated galaxies are at redshift 
z = 2. 



3 OBSERVATIONAL DIAGNOSTICS TO DISTINGUISH 
BETWEEN STAR FORMATION MODES 

In this section we present multiple observational diagnostics that 
can distinguish between quiescently star-forming galaxies (for 
SMGs this includes galaxy pairs and isolated discs) and starbursts 
induced at the final coalescence of major mergers. We present di- 
agnostics that rely only on integrated broadband photometry and 
CO line intensities (to determine gas mass). We also assume that 
sufficiently accurate redshifts are known. Spatially resolved data, 
such as that provided by (sub)mm interferometers and near-IR IFU 
spectrographs, can potentially provide more diagnostic power but 
come at a much greater observational cost. However, even with 
high-resolution near-IR data it can be difficult to distinguish be- 
tween disc galaxies and gas-rich mergers in which discs re-form 
shortly after final coalescence (Roberts on & Bullo ck 2008), and 
because of the high attenuation of SMGs near-IR observations may 
not probe the central starburst regions. The FIR and (sub)mm are 
able to probe much deeper into the central regions, but the spa- 
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Table 1. Simulation parameters 





M 2 oo a 


M*,init* 


]\/f . . c 
-^gas^nit 


-Rinit 


-ftperi 


el 


</>! 


el 


<\>1 


Name 


(lO 12 /^ 1 M ) 


(lO 10 ^- 1 M ) 


(lO 10 /*- 1 M ) 


(h' 1 kpc) 


(h' 1 kpc) 


(deg) 


(deg) 


(deg) 


(deg) 


b6i 


6.2 


5.3 


22 


70 


6.7 








71 


30 


b6j 


6.2 


5.3 


22 


70 


6.7 


-109 


90 


71 


90 


b6k 


6.2 


5.3 


22 


70 


6.7 


-109 


-30 


71 


-30 


b61 


6.2 


5.3 


22 


70 


6.7 


-109 


30 


180 





b6m 


6.2 


5.3 


22 


70 


6.7 








71 


90 


b6n 


6.2 


5.3 


22 


70 


6.7 


-109 


-30 


71 


30 


b6o 


6.2 


5.3 


22 


70 


6.7 


-109 


30 


71 


-30 


b6p 


6.2 


5.3 


22 


70 


6.7 


-109 


90 


180 





b5.5i 


3.2 


2.7 


11 


57 


5.3 








71 


30 


b5.5j 


3.2 


2.7 


11 


57 


5.3 


-109 


90 


71 


90 


b5.5k 


3.2 


2.7 


11 


57 


5.3 


-109 


-30 


71 


-30 


b5.51 


3.2 


2.7 


11 


57 


5.3 


-109 


30 


180 





b5.5m 


3.2 


2.7 


11 


57 


5.3 








71 


90 


b5.5n 


3.2 


2.7 


11 


57 


5.3 


-109 


-30 


71 


30 


b5.5o 


3.2 


2.7 


11 


57 


5.3 


-109 


30 


71 


-30 


b5.5p 


3.2 


2.7 


11 


57 


5.3 


-109 


90 


180 





b5i 


1.6 


1.4 


5.6 


44 


4.0 








71 


30 


b5j 


1.6 


1.4 


5.6 


44 


4.0 


-109 


90 


71 


90 


b5k 


1.6 


1.4 


5.6 


44 


4.0 


-109 


-30 


71 


-30 


b51 


1.6 


1.4 


5.6 


44 


4.0 


-109 


30 


180 





b5m 


1.6 


1.4 


5.6 


44 


4.0 








71 


90 


b5n 


1.6 


1.4 


5.6 


44 


4.0 


-109 


-30 


71 


30 


b5o 


1.6 


1.4 


5.6 


44 


4.0 


-109 


30 


71 


-30 


b5p 


1.6 


1.4 


5.6 


44 


4.0 


-109 


90 


180 






a Virial mass of each progenitor. 

b Initial stellar mass of each disc. 

c Initial gas mass of each disc. 

d Initial separation of the discs. 

e Pericentric passage distance. 

f Orientation of each disc's spin axis in spherical coordinates. 



tial resolution available at these wavelengths is much coarser un- 
less an interferometer is used, so typically only integrated FIR 
and (sub)mm photometry are available for high-redshift galaxies. 
(However, ALMA will soon change this situation drastically.) Fi- 
nally, mergers and starbursts are not equivalent (see above), so iden- 
tifying an object as a merger based on morphology does not en- 
sure it is also a starburst. Thus diagnostics that make use of only 
integrated data will continue to be crucial for distinguishing be- 
tween star formation modes and understanding the properties of 
high-redshift galaxies. 

We have identified the starburst phase by defining the baseline 
SFR as the minimum SFR that occurs between first passage and co- 
alescence and selecting all snapshots where the instantaneous SFR 
is > 3 x that baseline SFR. This factor is chosen so that the star for- 
mation induced by the merger dominates that which would occur 
in the discs even if they were not merging. Increasing (decreas- 
ing) the threshold would result in less (more) sources identified as 
bursts and amplify (diminish) the differences between modes that 
we describe below. The snapshots that meet the starburst criterion 
are labelled 'starburst' and plotted as blue squares. Since the mu- 
tual gravitational torques are sub-dominant at first passage relative 
to internal instabilities, the galaxies are primarily quiescently star- 
forming prior to the starburst induced at coalescence. We thus label 
all snapshots before the starburst phase 'quiescent' and plot them as 
black circles. All snapshots after the starburst phase are neglected, 
as these are typically AGN or spheroids with relatively little on- 



going star formation. We will investigate the FIR properties of ob- 
scured AGN in detail in future work. 

Our focus is the bright SMG population, defined by Ssso > 5 
mJy, so throughout this work simulated SMG data points are plot- 
ted with larger symbols (black circles for quiescently star-forming 
and blue squares for starburst SMGs). However, in order to com- 
ment on the applicability of the diagnostics to the high-redshift 
ULIRG population in general, data for simulated galaxies that 
would not be selected as SMGs are also plotted. Small black circles 
(blue squares) correspond to quiescently star-forming (starburst) 
ULIRGs with S 8 50 < 5 mJy. 

We have neglected all snapshots with greater than 40 per cent 
gas fraction in order to eliminate the early parts of the simulations 
when the discs are still stabilising and to remain consistent with (al- 
beit uncertain; e.g.,[B othwell et al. 2010) observational constraints 
(Tacconi et al. 2006 2008). Such high initial gas fractions are re- 
quired to maintain sufficient gas until the time of coalescence be- 
cause our simulations do not include any additional gas supply be- 
yond what the galaxies start with. We have checked that the results 
are qualitatively the same when we use an initial gas fraction of 
60 per cent and include all snapshots. All quantities plotted are 
totals for the entire system because we wish to present observa- 
tional diagnostics based on integrated data alone. In all figures, if a 
simulated galaxy sub-population has sample size N > 100 only a 
randomly selected subsample of 100 objects is plotted; this is done 
to improve the readability of the plots, as a sub-sample of 100 is 
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sufficient to show the distribution of the simulated data. Finally, 
note that we plot data from idealised simulations without applying 
any weighting to account for cosmological abundances. Thus the 
exact distribution of data in the various diagnostic plots is not nec- 
essarily representative of the SMG population. What is meaningful, 
however, is when sub-populations occupy distinct regimes in a di- 
agnostic plot; this is a clear prediction for how the star formation 
modes should differ and how one can observationally disentangle 
the classes in order to determine their relative contributions to a 
given galaxy population. 



3.1 Luminosity-effective dust temperature relation 

Far-IR (FIR) galaxy SEDs are often described in terms of a 'dust 
temperature' obtained via fitting a single-temperature (single-T) 
modified blackbody to the FIR SED (Hil debrand[T983] l. The single- 
T modified blackbody form is 



Su 



?o(l 



•e-<W)B v (T d ), 



(1) 



where S u is the flux density at rest-frame frequency v, So is the 
normalisation, Ai is the wavelength at which the effective optical 
depth r\ = 1, f3 is the effective slope of the emissivity in the FIR, 
Td is the effective dust temperatur^j and B u (Td) is the Planck 
function. The form is motivated by a simple model: If it is assumed 
that the entire mass of dust has physical temperature equal to Td, 
the source has constant source function B u (Td). Assuming also 
that the source is uniform, with projected area A and angular size 
dQ, and that the opacity curve is a power law in the FIR, 

Kx = k (A/Ao) _/3 , 



the optical depth is 



T\ = 



K \M d 

A 



KpM d ( A 
dQD 2 



(£)" 



(2) 



(3) 



where D is the distance to the source, Md is the dust mass, and ko 
is the opacity at wavelength Ao[j(For simplicity we have assumed 
z << 1 here. In practice, we transform our simulation photome- 
try to the rest frame before fitting.) Under these assumptions, the 
solution of the radiative transfer equation is 



S u = dQ(l- e 



dn 



1 — exp 



B v (T d ) 

-KoM d 
dQD 2 



(B 



B v {T d 



(4) 



(5) 



Note that the validity of Equation {3} is limited to the IR wave- 
lengths where the opacity curve has the assumed power-law form. 
Typically it is assumed that optical depths in the FIR are small, so 

[1 - exp (-(A/Ai)"")] » (A/Ai)"". Thus 



Su — S \ B v (Td), 



(6) 



6 Note, however, that in both the simulations and reality there is always 
a distribution of dust temperatures, as the temperature of a grain depends 
on the local radiation field and the size of the grain. The effective dust 
temperature fit in this manner may provide a good approximation to the 
luminosity -weighted average dust temperature, but to our knowledge this 
correspondence has not been demonstrated. 

7 It is important to distinguish between Ao and Ai . The former is just part 
of the parameterisation of the opacity curve of the dust (see Equation[2j. The 
latter depends on the dust mass and effective area of the source in addition 
to the dust opacity. Ai can be obtained from Equation |3j by setting t\ = 1 
and solving for A. 



where S = SoX^. We shall refer to the form given in Equation 
{6} as the single-T optically thin (OT) modified blackbody. We re- 
fer the reader to the appendix of Hll for more details about the 
modified blackbody forms. 

The simple model described above provides motivation for the 
fitting forms given in Equations {I} and {6}. However, because the 
simple assumptions inherent in the model (specifically the single 
dust temperature and assumed geometry) are certainly not true for 
real galaxies, one should not assume a priori that the parameters 
derived via fitting FIR SEDs with Equations |TJ or J6l can be used 
to infer the physical quantities present in Equation J5J. Indeed, the 
parameters determined using this fitting method should not be in- 
terpreted physically for multiple reasons: For example, (3 and Td 
are degenerate (e.g., Sajina et al. 2006), and the Td and f3 one 
derives from the fitting depends strongly on observational noise 
(Shet ty et al.||2009a|), dust temperature variations along the line- 
of- sight (Shetty et al. 2009b), and the wavelength range spanned 
by the photometry included in the fit ([Magnelli e t al. 12012] here- 
after Ml 2). Even in the case of a uniform radiation field grains 
of different sizes will have different temperatures, so there is no 
reason to expect (3 derived from the single-T modified blackbody 
to be equal to the intrinsic power-law index of the dust emissiv- 
ity in the FIR. (However, if a distribution of dust temperatures is 
permitted it is perhaps possible to recover the true f3 from fitting 
the SED.) Furthermore, adding a significant component of very 
cold dust can mimic the effect of high effective optical depth, so r 
and the temperature distribution are also degenerate (Papadopoulos 
|et al.|2010b] ). Finally, even the fitting method can cause incorrect 
conclusions: when traditional non-hierarchical % 2 minimisation is 
used, observational noise can introduce a spurious anti-correlation 
between Td and f3 even when the true values of Td and f3 are posi- 
tively correlated; it is possible to avoid this pitfall by utilising a hi- 
erarchical Bayesian method ( |Kelly et al.|2012| ). Thus, at the least, 
one should use caution when attempting to infer physical condi- 
tions from the parameters derived by modified blackbody fitting, as 
we demonstrate below. 

However, even if the model parameters cannot be interpreted 
physically, the modified blackbody forms often provide acceptable 
descriptions of the FIR SEDs of galaxies. Thus in order to compare 
to observations in a meaningful way and provide testable predic- 
tions from our models we have fit the FIR SEDs of our simulated 
galaxies with both the optically thin and full forms of the modified 
blackbody. We have also used a more sophisticated model which 
assumes a power-law distribution of temperatures. The results for 
each fitting form are discussed in turn below. 



3.1.1 Single-T, OT modified blackbody 

Fig. [3] shows the distribution of our simulations on the Td — Tir 
plot when the single-T OT modified blackbody form (Equation [6} 
is used. The top panel shows Td derived from fitting the simu- 
lated Her schel SPIRE (Griffin et al.|2010) 250-, 3 50-, and 500- urn, 
SCUB A ( [Holland et allT999] ) 850-^im, and AzTEC ( [Wilson et al.| 
2008) 1.1 -mm photometry; on the bottom the simulated Herschel 
PACS (Poglitsch et al. 2010) 160-pjn photometry are also included 
in the fit. We have excluded the PACS 70- and lOO-i^m points be- 
cause often an acceptable fit cannot be obtained if those points 
are used. Possible reasons for this are that, at those wavelengths, 
stochastically heated grains contribute significantly to the SED and 
the assumption of optical thinness is violated most. When perform- 
ing the Levenberg-Marquardt least-squares fit we have assumed 10 
K < T d < 100 K, 1.0 < P < 2.5, and constant fractional flux er- 
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Table 2. Median T d for the OT single-T modified blackbody 



Galaxy type 



without PACS 160-um 

L m < 10 13 L L m ^ 10 13 L 



with PACS 160-um 

L IR < 10 13 L L m ^ 10 13 L 



sample size 

L IR < 10 13 L L IR ^ 10 13 L 



Quiescent SMGs 


30.1 K 


36.4 


34.2 


42.0 


954 


8 


Starburst SMGs 


33.4 


38.1 


35.1 


41.0 


409 


292 


Quiescent non-SMGs 


28.9 


- 


34.1 


- 


648 


- 


Starburst non-SMGs 


35.2 


42.0 


36.4 


50.2 


741 


21 



70 



> quiescent SMGs 
i starburst SMGs 

quiescent non-SMGs 
■ starburst non-SMGs 



* 50 




70 p 



o quiescent SMGs 
□ starburst SMGs 
• quiescent non-SMGs 
- starburst non-SMGs 



with PACS 



* 50 = 



40; 



30 = 




10 13 
L, R (U) 



10 14 



Figure 3. Effective dust temperature Td (K) derived from fitting the single- 
T OT modified blackbody (Equation [6} to the simulated SPIRE, SCUBA, 
and AzTEC photometry (top) and PACS 160-um, SPIRE, SCUBA, and 
AzTEC photometry (bottom) versus total IR luminosity Lir ( L ). For 
the SMG population, effective dust temperature correlates strongly with lu- 
minosity, and the sources in the high-Td, high-LiR region of the plot are 
almost exclusively merger-induced starbursts. The SMG selection is biased 
against galaxies with hotter SEDs, so, for a given luminosity, the non-SMGs 
have higher T& than the SMGs. At the lower luminosities there is a signif- 
icant population of galaxies missed by the SMG selection; these cause the 
Td — Lm relation for the full population to have larger scatter than that for 
the SMG population. Note also that adding the PACS 160-um data results 
in higher T^ and increased scatter. 



(35) K are starbursts. Thus one can use a cut in Td to cleanly se- 
lect starbursts from the overall SMG population. Furthermore, for 
the SMG population the cut Lir > 10 13 L can also be used to 
cleanly select starburst SMGs. 

When the non-SMG population is also considered the scatter 
in the Td — Tir relation is increased, as there is a significant popu- 
lation of galaxies at the lower-LiR end with relatively high Td that 
are missed by the SMG selection (see also fig. 1 of M12). These 
simulated galaxies correspond to the observed 'hot-dust ULIRGs' 
([Chapman et al.||2004l [Ch apman et al. 2008] |Casey et al.|[2009| 
|Casey et al.|2011[|Magnelli et al.|20To[ M12); we discuss the dif- 
ferences between the SMG and non-SMG sub-populations in detail 
below (Section [43) . At a given Lir the non-SMGs have higher 
effective dust temperatures than the SMGs because of the bias in- 
herent in the SMG selection[jAs for the SMG population, almost 
all non-SMGs with T d > 40 (35) K when the PACS 160-um point 
is (not) used in the fit are starbursts. Thus the Td cut is an effective 
means for selecting starbursts from the overall galaxy population 
even at luminosities where there are a mix of starbursts and qui- 
escently star-forming galaxies. This is consistent with the observa- 
tions of |Magnelli etaTlpOTO] ) and M12. 

Note also that inclusion of the PACS photometry results in 
both increased dust temperature and larger scatter. This effect oc- 
curs because we allow f3 to vary; if f3 is fixed then Td is relatively 
similar whether or not the PACS photometry is included (Magdis 
|et al.|2010[ M12). The fact that the results of IR SED fitting can 
depend strongly on the regions of the SED sampled is another rea- 
son one must use caution when interpreting the derived parame- 
ters. This is especially significant when comparing low- and high- 
redshift samples, as in this case a given observed data point can 
correspond to significantly different rest-frame wavelengths for dif- 
ferent redshifts. 

In Fig. [4] we compare our simulations to the observations of 
Ml 2. The effective temperature plotted here has been determined 
by fitting the single-T OT form with fixed /3 = 1.5, which is 
what |Magnelli et al.| have done. Only the PACS 160-um, SPIRE, 
SCUBA, and AzTEC photometry have been used for the fits. For 
the range in Lir covered by the simulations the agreement between 
the simulated and observed SMGs is very good. This agreement 
gives confidence that the FIR SEDs of the simulated galaxies are 
reasonable. However, there are a few galaxies that lie below the re- 
gion spanned by the simulations. This is consistent with the conclu- 
sion of Jonsson et al. (2010) that real galaxies form a more diverse 
family of SEDs than are generated by the simulations. 



ror of ten per cent. The median effective dust temperatures for the 
each galaxy type are given in Table[2] where the galaxies have been 
divided into two bins, Lir < 10 13 L and Lir ^ 10 13 L . 

The key trends to take away from Fig. [3] are that effective 
dust temperature correlates with luminosity, and the most lumi- 
nous, hottest SMGs are almost exclusively starbursts. When the 
PACS photometry are (not) used, almost all SMGs with T d > 40 



8 Note that this is not necessarily reflected in the tables of median temper- 
atures because the non-SMGs in a given bin tend to have lower Lir than 
the SMGs; thus the comparison of the median Td values for the bins is not 
the same as comparing the values for populations at fixed Lir. 
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> quiescent SMGs 
i starburst SMGs 

quiescent non-SMGs 
■ starburst non-SMGs 

> M12 unlensed 

> M12 lensed 



10 11 




> quiescent SMGs 
i starburst SMGs 
quiescent non-SMGs 
starburst non-SMGs 



40E-\tfj 



without PACS 



L, R (U) 




Figure 4. Effective dust temperature (K) yielded by fitting the single-T OT 
form of the modified blackbody with /3 = 1.5 fixed versus Lir ( L©). In 
order to compare directly to Ml 2 only the PACS 160-jxm, SPIRE, SCUBA, 
and AzTEC photometry have been used in the fit. The large black circles 
(blue squares) are the simulated quiescently star-forming (starburst) SMGs, 
and the small black circles (blue squares) are the simulated quiescently star- 
forming (starburst) non-SMGs. The red diamonds (green triangles) with 
error bars are the unlensed (lensed) SMGs (Sg50 > 3 mJy) from M12. 
The agreement between the simulated and observed unlensed SMGs is very 
good for the entire luminosity range sampled by our simulations, giving 
confidence that the FIR SEDs of the simulated galaxies are reasonable. 



3.1.2 Full form of the single-T modified blackbody 

There is increasing evidence that the simple single-T OT modified 
blackbody form given in Equation (6) provides a poor fit to the FIR 
SEDs of simulated (Hll) and observed ( |Lupu et a l. 2010, Conley 
[eTaL][20TT| M12; Sajina et al., submitted) high-redshift ULIRGs 
on the Wien side of the SED. Instead, more sophisticated forms, 
such as Equation {T} or multiple-temperature models (e.g., Dale & 



|Helou|2002l|Clements, Dunne, & Eales|2010| [Kovacs et al.|20To^ 
must be used. Equation (I} allows non-negligible optical depths in 
the FIR, and allowing f3 to vary mimics the effect of a tempera- 
ture distribution (e.g., |Shetty et al.|2009b[|Clements et al.|2010) . 
The multiple-temperature models can account for non-negligible 
optical depths in the IR and multiple temperatures of dust, both of 
which are physically more valid assumptions than those implicit in 
the single-T OT blackbody model, but require at least one parame- 
ter more than Equation ([TJ. 

Fig. [5] shows the Td — Tir plot when we derive the effec- 
tive dust temperature Td by fitting the full form of the modified 
blackbody, Equation ([TJ, to the simulated SPIRE+SCUBA+ AzTEC 
(top) and PACS+SPIRE+SCUBA+AzTEC (bottom) photometry. 
The PACS 70-pjn data have not been used because at z — 2 
observed-frame 70 |^m corresponds to rest-frame 23 \xm, a re- 
gion of the SED dominated by emission from stochastically heated 
grains. If the 70- |jjn data are included an acceptable fit is often not 
possible. We have assumed 10 K < T d < 100 K, 1.0 < /3 < 2.5, 
and ten per cent flux uncertainty. Ai is allowed to vary freely, but in 
practice it is always greater than ^375 ^im. The median effective 
dust temperatures are given in Table [3] 

Almost all trends demonstrated by Fig.[3]hold here: effective 
dust temperature correlates with luminosity, and the most lumi- 
nous, hottest sources (Td > 70 (60) K when the PACS data are 
(not) used) are almost exclusively starbursts. Almost all galaxies 
with Lir > 10 13 L are starbursts, but, again, a Td cut is sig- 
nificantly better than an Lir cut at selecting starbursts from the 



()() 








1 o quiescent SMGs 


with PACS E 




: a starburst SMGs 


- 


90 


■ quiescent non-SMGs 


: 




; ■ starburst non-SMGs 


□ E 


80 


r a ft n a J 1 □ 


D _E 




: „" " n D Oft 


E 






: 


/O 




- 




: ■ ■ „ ■ o^bfl n no D 


= 


60 
50 


-E 




m n 


- 


40 


-~ 


-E 




, , , , 1 


, , , , - 



10 1 ' 



10 13 
L, R (U) 



10 14 



Figure 5. Same as Fig.|3] but the effective dust temperature has been derived 
by fitting the full form of the modified blackbody (Equation [TJ. The Td 
values derived using this form are systematically higher than when optical 
thinness is assumed (see text for details), but all trends - except for the 
increased scatter when PACS 100- and 160-um data are included - seen in 
Fig.[3]still hold. 

non-SMG population. There is a significant population of galaxies 
at the lower luminosity end that are missed by the SMG selection 
because of their relatively high effective dust temperatures. Finally, 
inclusion of the PACS photometry again results in higher effective 
dust temperature. 

Comparison of Figs. [3] and [5] and the median values given 
above shows that assuming optical thinness results in systemati- 
cally lower Td than when Equation {I} is used. This occurs be- 
cause [l - exp (-(A/Ai) -/3 )] < (A/Ai) -/3 for all A > 0. Thus, 
for fixed Td and f3, the assumption of optical thinness will system- 
atically over-predict the flux at frequencies for which r > 1. As 
a result, for a given SED and fixed f3, Td derived from Equation 
(61 will be lower than that derived from Equation (pi. This effect 
has been demonstrated when fitting SEDs of high-redshift ULIRGs 
( [Lupu et al.|2010||Conley et al.|20lT| Sajina et al., submitted), and 
it shows that one should use caution when attempting to interpret 
Td physically. 

3.1.3 Power- law T -distribution model 

We have also fit a subset of the simulated galaxies' SEDs assuming 
a power-law temperature distribution with a low-temperature cutoff 
following Kovacs et al. (2010). The fitting method is summarised 
here, but the reader is referred to Kovacs et al. (2010) and Ml 2 for 
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Table 3. Median T& for the full form of the single-T modified blackbody 



Galaxy type 



without PACS 100- & 160- urn 

L m < 10 13 L L m ^ 10 13 L 



with PACS 100- & 160-M-m 
L IR < 10 13 L L IR ^ 10 13 L 



sample size 

L IR < 10 13 L L IR ^ 10 13 L 



Quiescent SMGs 


42.7 K 


63.5 


60.4 


69.6 


954 


8 


Starburst SMGs 


48.6 


61.3 


57.2 


70.4 


409 


292 


Quiescent non-SMGs 


40.9 


- 


60.1 


- 


648 


- 


Starburst non-SMGs 


53.3 


80.1 


60.0 


77.8 


741 


21 



Table 4. Median low-temperature cut-off T c for the power-law dust temperature distribution model 





with PACS 100 


- & 160-^tm 






sample size a 




Galaxy type 


L m < 10 13 


L 


L m > 10 13 


L 


Lir < 10 13 L 


L m ^ 10 13 L 


Quiescent SMGs 


31.0K 




34.7 




31 






1 


Starburst SMGs 


32.6 




40.9 




41 






46 


Quiescent non-SMGs 


27.6 




- 




15 






- 


Starburst non-SMGs 


36.0 




- 




164 






- 



2 The sample sizes here differ from those above because only a subset of the simulated galaxy SEDs have been 
fit with the power-law T- distribution form. 
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Figure 6. Same as Fig. [3] but the effective dust temperature has been de- 
rived by assuming a power-law distribution of dust temperatures following 
Kovacs et al. (2010, the temperature plotted is the low-temperature cutoff 
of the power-law distribution - see text for details). Again, the dust tem- 
peratures differ, but the trends are insensitive to the manner in which the 
effective dust temperature is derived. All sources with T c > 36 K are star- 
bursts. 



full details of the model. In this model, the dust has a distribution 
of physical temperatures given by dMd/dT oc T -7 for T > T c 
and Md(T) = otherwise. The effective optical depth t\ is pa- 
rameterized as 



T\ = Kq 



(*)' 



M d 



(7) 



where R e ^ can be thought of as an effective radius of the source. 
Because of the added parameter 7 we have always used the PACS 
100- and 160-^im, SPIRE, SCUBA, and AzTEC data. Following 
M12, we have assumed that single values of f3, 7, and i? e ff can be 
used for all sources, and we have used a subset of 20 simulated 
galaxies to fix those parameters in the following manner: We first 
gridded the (/?, 7, R G s) parameter space. For each point in the grid 
we fit all 20 sources allowing Td and Md to vary. We summed the 
X 2 values of the individual fits for each parameter combination and 



chose the parameter combination with the lowest total % 2 value. As 
above, the fractional flux error assumed is ten per cent. The param- 
eters we determined in this manner are (/3, R G s , 7) = (1.6, 2 kpc, 
8.7). |Kovacs et al. (2010) found that the parameters (/?, R e & , 7) = 
(1.5, 1 kpc, 6.7) provided the best description of their sample of 
z ~ 2 starbursts. For the SMG sample of M12 the best-fitting pa- 
rameters are (/3,# e ff,7) = (2.0 ± 0.2, 3 ± 1 kpc, 7.3 ± 0.3). 
The values of ft and i? e ff derived from our simulations lie between 
those from the two observational studies. However, our temperature 
distribution is steeper than those found by both Kovac s~e"t al.| and 
Magnelli et al. perhaps because our simulations do not yet include 
stochastically heated very small grains and thus may underestimate 
the amount of dust at high temperatures. We defer a detailed com- 
parison of simulated and observed ULIRG SEDs and the derived 
dust parameters to future work. 

The resulting T c — Lm plot is shown in Fig. [6] Note that the 
temperature plotted here is the low-temperature cutoff, which is 
also the temperature of most of the dust because of the steepness of 
the power-law distribution. The median T c values are given in Table 
HI As seen for both single-T fitting forms, for SMGs there is a clear 
correlation between effective dust temperature and luminosity, and 
the starbursts are the most luminous and have the highest values 
of T c . An effective dust temperature cut is very effective here: all 
galaxies with T c > 36 K are starbursts. Furthermore, at the lower- 
luminosity end there is again a population of starbursts missed by 
the SMG selection because of their relatively high effective dust 
temperatures. Consequently, a T c cut can select a larger subset of 
the starbursts than can an Lir cut. 

M12 argue that the cut T c > 25 K can very effectively sep- 
arate starbursts (identified observationally by their offset from the 
SFR-M* relation) from quiescently star-forming SMGs, as we also 
see in our simulations. However, the value of T c that separates sim- 
ulated starburst SMGs from quiescently star-forming SMGs differs 
significantly from that derived by Ml 2. The primary reason for this 
discrepancy is that the (/5, i2 e ff , 7) values derived for our simula- 
tions differ from those for the M12 sample. If the simulated SEDs 
are fit with the (/3, i? e ff , 7) values from Ml 2, the resulting T c val- 
ues are significantly lower, the SED fits are still acceptable, and 
the agreement between the observations and simulations is much 
improved. See section 6 of Ml 2 for further details. 
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Table 5. T d and (3 for the best-fitting models 
shown in Fig.pl 



Model 



P 



Single-T OT form, fixed f3 47 K 1.5 

Single-T OT form 53 1.1 

Single-T full form 69 1.5 

Power-law T-distribution 49 1.2 
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Figure 7. Effective dust temperature derived from fitting the single-T OT 
modified blackbody (Equation [6] top), full modified blackbody (Equa- 
tion [^ middle), and power-law T-distribution model (bottom) to the 
PACS+SPIRE+SCUBA+AzTEC integrated photometry versus nuclear sep- 
aration for snapshots where the BHs have not coalesced. T d is anti- 
correlated with nuclear separation, but there is a large scatter in effective 
temperature at a given nuclear separation. The objects with the highest ef- 
fective dust temperatures (T d;0 T > 35 K, T d > 55 K, and T c > 36 K) 
are almost exclusively late- stage merger-induced starbursts with 4h ^5 10 
kpc. The simulations with high T d at 4h ~ 100 kpc are those where Lir 
- and thus T d - is still high (but instantaneous SFR is not) because of the 
burst of star formation that can occur at first passage. 



3.1.4 Summary of results that are independent of the fitting form 

Figs. [3] [5] and [6] all show that the T d — Tir plot is an excellent 
way to select starburst SMGs from the general population. In all 
three figures there is a clear correlation between Td and Tir, which 
agrees with observations of both local and high-redshift ULIRGs 
(e.g.,|Kovacs et al.|2006[|M agnelli e tal.|2010[ Amblard eTal.|2010[ 
IChapman et al.|2010[|Hwang et al.|2010[ lMagdis e tal.|2010[ M12). 



Though there is some overlap between the quiescently star-forming 
and starburst sub-populations, the most luminous, hottest SMGs 
are almost exclusively starbursts. Note that both the T d — Tir cor- 
relation and the separation between the populations are indepen- 
dent of the fitting method used, though the specific temperature 
values above which there are no quiescently star-forming galaxies 
differ (as expected because of the systematic difference in temper- 
atures yielded by the two methods). Thus our simulations make the 
clear, robust prediction that the most luminous galaxies will have 
the hottest SEDs and will typically be late-stage merger-induced 
starbursts. 

In all cases there is a sub-population of hot-dust ULIRGs 
which have relatively high T d for a given IR luminosity. Such 
galaxies are not present in the SMG sub-population because of the 
bias of the SMG selection. The SMG selection bias results in an 
apparent relatively tight correlation between T d and Tir, so an 
Tir cut is roughly as effective as a T d cut at selecting a subset of 
starburst SMGs from the general SMG population. When the hot- 
dust ULIRGs are included the scatter in the T d — Tir relation for 
Tir < 10 13 L is increased significantly. For Tir < 10 13 L a 
T d cut can select starbursts whereas an Tir cut cannot. Thus if one 
wishes to select starbursts from a given galaxy population a T d cut 
is preferred because it can isolate a larger subsample of starbursts 
with a wider range in Tir. 

Fig.|7]shows T d derived from fitting the single-T OT modified 
blackbody form (top), full form of the single-T modified black- 
body (middle), and power-law T-distribution model (bottom) to the 
PACS, SPIRE, SCUBA, and AzTEC photometry (with some PACS 
data excluded for each form, as described above) versus separa- 
tion of the central BHs (aka nuclear separation; 4h). The nuclear 
separation serves as a proxy for the merger stage, but it is impor- 
tant to keep in mind that it does not decrease monotonically as 
the merger progresses. The starburst galaxies have systematically 
lower 4h than the quiescently star-forming galaxies because the 
tidal torques which drive the starburst are strongest at final coales- 
cence of the two discs. The typical T d values increase with decreas- 
ing c/bh, though there is large scatter at a given dsH, especially for 
4h ^$10 kpc. This occurs because, for a given simulation, Tir is 
anti-correlated with 4h (because the starburst is strongest at low 
nuclear separations), and, as seen above, the most luminous galax- 
ies are also the hottest. Interestingly, for a given nuclear separation 
the SMGs and non-SMGs span a similar range in T d . The simu- 
lated mergers with 4h ^100 kpc and relatively high T d are those 
observed shortly after first passage. At this time Tir is high and 
the starburst mode is important even though the galaxies do not 
meet the definition of starburst because the instantaneous SFR has 
dropped. 
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Figure 8. Top: The black line is an example rest-frame FIR SED, \L\ 
( L0 ) versus A res t (M-m), of a simulated starburst. The crosses are the PACS 
100- and 160- urn, SCUBA 250-, 350-, and 500-um, SCUBA 850-um, and 
AzTEC-1.1 mm photometry calculated from the simulation SED. The other 
lines are best fits to the photometric points for different fitting forms: the 
single-T, OT modified blackbody (Equation [6j with fixed (3 = 1.5 (ma- 
genta dashed), the single-T OT modified blackbody with (3 allowed to vary 
(red long-dashed), the full form of the single-T modified blackbody (Equa- 
tion[T] blue dash dot), and the power-law T-distribution ( Kovacs et al. 2010 , 
green dash dot dot dot). Bottom: Ratio of the fitted L \ to the actual L \ ver- 
sus A r est (pn). The crosses overlaid on the black solid line, which denotes 
L\,fit = L\, show the wavelengths of the photometric points used for the 
fits. The values of T^ and (3 for the best- fitting models are given in Table|5] 
All forms except the single-T OT form with fixed (3 recover the photome- 
try within the 10 per cent uncertainty. However, the derived dust tempera- 
ture and (3 depend strongly on the fitting method, so the model parameters 
should not be interpreted physically even though the models adequately de- 
scribe the data. 



when performing the \ 2 minimisation. The single-T OT form with 
fixed /3 can only recover the photometry to within ~ 30 per cent. 
Though the more complicated forms are successful at recovering 
the photometry used for the fitting, they have varying levels of suc- 
cess describing the SED beyond the wavelengths spanned by the 
photometry. As explained above, we expect the OT modified black- 
body to under-predict the SED on the Wien side of the SED. In- 
deed, this model under-predicts the SED shortward of the lOO-pm 
point, and the under-prediction is more severe when /3 is fixed. The 
full form of the modified blackbody fares better, but it also under- 
predicts the SED for rest-frame wavelength A res t < 20 |^m. The 
power-law temperature distribution model fares best at the shortest 
wavelengths, but it over-predicts the SED at A res t ~ 15 — 25 \im. 
At the longest wavelengths the single-T OT and power-law T- 
distribution models are most accurate because they have relatively 
low (3 values and thus less steeply declining SEDs. 

Table [5] shows that the derived parameters Td and f3, which 
are often interpreted as a physical dust temperature and the intrin- 
sic power-law index of the emissivity of the dust grains in the FIR, 
vary for the different fitting forms. Since observed galaxies typi- 
cally have 20 < T d < 80 K and 1 < /3 < 2, the variation among 
the best-fitting values is very significant. Consequently, it is diffi- 
cult to interpret the fitted parameters physically, as clearly the in- 
trinsic properties of the dust do not vary with the method used to fit 
the SED [j If the fitted effective dust temperatures have a physical 
meaning, they may correspond to different physical temperatures 
(e.g., the single-T modified blackbody may recover the luminosity- 
weighted dust temperature whereas the power-law model may bet- 
ter recover the mass- weighted temperature). None of the fitted f3 
values recover the intrinsic /3 of the dust, which is ~ 2 for the dust 
model we use. This is not unexpected, because a distribution of 
physical dust temperatures will change the slope of the SED and 
thus the f3 in a single-T model; non-negligible optical depths in the 
IR further complicate the picture. 

Clearly it is necessary to determine how the fitted parameters 
relate to intrinsic properties of the dust, but we defer further explo- 
ration of this complex topic to future work. We only wish to stress 
that it may be necessary to use forms more sophisticated than the 
single-T OT modified blackbody to fit IR SEDs and that the pa- 
rameters derived from the fits should not be interpreted physically. 
Instead, the models should be thought of as useful ways to encapsu- 
late the data with a few parameters and thus compare galaxy SEDs 
in a simple way by comparing the fitted parameters. Put another 
way, the Td — Tir plot still contains useful information about SED 
variation even if Td is not a physical dust temperature, and differ- 
ences in Td amongst galaxies reflect real differences in the galax- 
ies' SEDs. 



3.1.5 Comparison of fitting forms 

In the top panel of Fig. [8] we show the rest-frame SED of one of 
the simulated starbursts viewed from a single viewing angle. The 
over-plotted data points are the simulated photometry as described 
in the caption. The four lines are fits to the photometry using the 
models discussed above. The bottom panel shows the ratio of the 
model SED derived from fitting the photometry to the actual SED 
for each of the fitting methods. The values of Td and f3 for the best- 
fitting models are given in Table [5] 

All forms except the single-T OT modified blackbody with 
fixed p — 1.5 are able to recover the simulated photometry to 
within ~ 10 per cent, which is the level of uncertainty assumed 



3.2 Star formation efficiency 

Since starbursts form stars much more efficiently than quiescently 
star-forming galaxies, one should be able to distinguish between 
them via some observational indicator of star formation efficiency. 
To be consistent with the literature (e.g., Daddi et al. 2010, Gen- 
|zel et al.||2010| ) we define 'star formation efficiency' as SFE = 
^iR/Mg as . Note, however, that Lir does not necessarily track the 



9 Note that the issues discussed here are independent of observational 
noise, as we have added none to our simulated SEDs. Observational noise 
further complicates interpretation of the derived parameters (Shetty ~et al.| 
[2009a][KeTr7et al.|2012| . 
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Figure 9. 'Star formation efficiency' SFE = LiR/M gas ( Lq/ Mq) versus 
Lir ( Lq; top) and nuclear separation (kpc; bottom). SFE and Lir are 
totals for the system in order to mimic the effects of blending in the FIR and 
(sub)mm. At fixed Lir the SFE is characteristically higher by a factor of 
a few for starburst SMGs than for quiescently star-forming SMGs, and the 
starbursts (quiescently star-forming galaxies) missed by the SMG selection 
have SFE greater than the starburst (quiescently star-forming) SMGs. The 
spread in SFE at fixed Lir between the most extreme starburst non-SMGs 
and the quiescently star-forming SMGs can be as much as ~ 10 x. The 
starbursts near coalescence have the highest SFE. The objects with high 
SFE at 4h ~ 100 kpc are the same as the high-Td objects with 4h ~ 
100 kpc, mergers observed shortly after first passage when Lir is still high 
but instantaneous SFR is not. 



Table 6. Median SFE values ( L / M ) 



Galaxy type 



Lir < 10 13 L Lir 3* 10 ld L G 



Quiescent SMGs 


13.6 


60.3 


Starburst SMGs 


37.1 


91.2 


Quiescent non-SMGs 


9.4 


- 


Starburst non-SMGs 


51.2 


123.4 



instantaneous SFR because, in addition to recently formed stars, 
older stars and AGN can also contribute to Lir (see, e.g., section 
2.5 of Kennicutt 1998a and section 3.1 of Hll for details), so this 
should be considered only an approximate indicator of star forma- 
tion efficiency. Furthermore, inferring the total molecular gas mass 
from CO observations is notoriously difficult, as the CO-H2 con- 
version factor Xco = Nu2/ ' L' co depends on the giant molecu- 
lar cloud surface density and the kinetic temperature and velocity 



dispersion within clouds (Narayanan et al. 2011b c, S hetty et al.| 
|2011a|b| |Papadopoulos et al.|2012| ). As a result, Xco is expected 
to be a factor of ~ 2 — 10 x lower in starbursts than in disc galaxies 
(Narayanan et al. 2011b). The uncertainty surrounding Xco com- 
plicates efforts to determine how much the SFE of starbursts and 
quiescently star-forming galaxies differs (e.g., Papadop oulos et al.] 
|2012| ). It would be best for us to predict the CO line luminosity for 
our simulated galaxies, as has been done by, e.g., Narayanan et al.| 
(2009), and, ideally, to self-consistently track formation and de- 
struction of molecular gas (see, e.g., Robertson & Kravtsov 2008). 
However, doing so requires introduction of another code - and the 
corresponding complexities and uncertainties - in addition to the 
two employed, so we feel this is best left to future work (though we 
briefly discuss how the molecular gas emission might be used as a 
diagnostic in Section |4~4] >. Thus when calculating SFE here we use 
the total gas mass instead of the molecular gas mass. 

In Fig. [9] we plot SFE versus Lir (top) and nuclear separation 
(bottom); the median values are given in Table [6] At a given Lir, 
the starburst SMGs have SFE up to 5 x greater than that of the qui- 
escently star-forming galaxies. The starbursts not selected as SMGs 
have the highest SFE, whereas the quiescently star-forming SMGs 
have the lowest values; the discrepancy can be as much as ~ 10 x . 
The bottom panel of Fig. [9] demonstrates that for both the SMGs 
and non-SMGs the SFE increases as the merger advances and is 
highest for mergers nearest coalescence. At a given 4h the SMGs 
and non-SMGs span a similar range in SFE. The objects with high 
SFE at 4h ^100 kpc are mergers observed shortly after first 
passage. For these galaxies Lir is still high, causing high SFE, but 
instantaneous SFR is not, so they do not meet the definition of star- 
burst. 

The simulations qualitatively reproduce the systematic offset 
in SFE between quiescently star-forming galaxies and starbursts 
shown in fig. 1 of |Daddi et al.| ( |20101 ). However, for a given Lir 
the SFE values of the simulated starbursts are significantly (~ 5 x ) 
lower than the observed values. Part of this is because our measure 
of the SFE in the simulations uses the total gas mass, not the molec- 
ular gas mass. In the simulations at the peak of the starburst, cold 
gas within 5 kpc of the centre, which we take as a rough approxi- 
mation for molecular gas that would be probed by observations (see 
Narayanan et al. 2009 ), is typically less than half of the total gas 
mass even though it accounts for effectively all of the star forma- 
tion. Thus if we were to calculate SFE for our simulated starbursts 
using molecular gas mass rather than total gas mass the values for 
the simulations would be a factor of ~ 2 — 3 x higher, which would 
account for a large part of the discrepancy between the simulated 
and observed starbursts. 

The SFE of the simulated starburst SMGs is only a factor 
of ~ 2 — 3 x greater than that of the simulated quiescently star- 
forming galaxies, whereas the observed difference is ~ 4 — 10 x. 
For the starbursts nearest coalescence, however, the difference can 
be as great as 10 x (see bottom panel of Fig. [9}. One possible rea- 
son the SFE discrepancy is lower than observed is that, because of 
the setup of our simulations, the pre-burst, quiescently star-forming 
discs have systematically higher gas fractions than those discs near 
coalescence. Since the star formation law implemented in the sim- 
ulations has a non-linear dependence on gas density the SFE will 
increase with gas fraction. Despite the possible discrepancies de- 
scribed above, the simulations make the robust prediction that those 
systems with largest SFE at a given Lir will be merger-induced 
starbursts. 
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Figure 10. IR excess (IRX = Lir/Lfuv) versus Lir ( L ; top) and nu- 
clear separation (kpc; bottom). Both IRX and Lir are calculated for the to- 
tal system. For the starbursts, IRX is correlated with Lir and anti-correlated 
with nuclear separation. The starbursts are typically more obscured (have 
higher IRX) than the quiescently star-forming galaxies. However, for fixed 
Lir the starbursts and quiescently star-forming discs have similar IRX val- 
ues, so the reason starbursts typically have higher attenuation is because 
they are typically more luminous than the quiescently star-forming discs. 
Interestingly, the correlation is reasonably well-approximated by constant 
Lfuv = 4 x 10 10 L© (solid line in the top panel) at the bright end. 



Table 7. Median IRX values 



Galaxy type 



Lir < 10 13 L Lir > 10 13 L 



Quiescent SMGs 


144.9 


340.0 


Starburst SMGs 


185.4 


332.3 


Quiescent non-SMGs 


120.1 


- 


Starburst non-SMGs 


199.9 


373.3 



3.3 IR excess 

In Fig. [To] we plot the total IR luminosity divided by the rest-frame 
far-UV luminosity, Lir/Lfuv - referred to as the IR excess (IRX) 
- versus Lir (top) and nuclear separation (bottom); the median val- 
ues are given in Table|7] Both IRX and Lir are calculated for the to- 
tal system. The IRX serves as a measure of the level of obscuration 
of a galaxy. IRX increases with Lir, as has previously been both 



observed (e.g., Wang & Heckman 1996, Buat & Burgarella 1998, 
Buat et al.|1999||2005||2b07. 2009, Adelberger & Steidel 2000, 



Figure 11. SFR ( M yr _1 ) versus M* ( M ), where both are totals for 
the system. The solid line is the relation for the z = 2.0 — 2.5 bin of Karim 
|et al.| (2011) . Simulated galaxies above the observed relation are almost 
exclusively starbursts, and the starbursts typically lie above the relation. 



Hopkins et al. 2001, Bell 2003, Reddy et al. 2010) and demon- 
x H 1 1 H 1 1 j ii i 



strated by simulations ( Jonsson et al. 2006). The quiescently star- 
forming galaxies tend to have lower IRX than the starbursts, pri- 
marily because the starbursts are typically more luminous. At a 
given Lir, the starbursts and quiescently star-forming discs have 
very similar IRX. Jonss o~et al.| d2006 ) have previously demon- 
strated this result using similar simulations, and they argue that the 
correlation arises because both SFR and dust optical depth correlate 
with density. The bottom panel shows that IRX increases as nuclear 
separation decreases. This adds further evidence in support of the 
arguments of Jonsson et al. (2006): The galaxies are more compact 
at coalescence than during the pre-burst, infalling-disc stage. The 
resulting higher densities lead to higher SFR and thus greater Lir. 
Furthermore, the stars formed in the starburst, which dominate the 
luminosity, are centrally concentrated and thus typically more ob- 
scured than stars distributed throughout the initial discs. Thus Lir 
and IRX both increase as nuclear separation decreases. This find- 
ing is consistent with that from studies of AGN obscuration using 
similar simulations ( [Hopkins et al.|2005a|b] 2006}. 

The solid line in the top panel of Fig.[l0|is the relation given 
by constant Lfuv = 4 x 10 10 L©. Note that, assuming this value, 
Lfuv is less than one per cent of the bolometric luminosity. At the 
bright end this relation approximates that of the simulated galaxies 
reasonably well, suggesting that the light observed in the rest-frame 
FUV is decoupled from the bolometric luminosity. This what we 
expect when the luminosity of a galaxy is dominated by a deeply 
dust-enshrouded starburst and the only light observed in the FUV is 
from stars located outside the heavily obscured nuclear region (see 
also Jonsson et al. 2006). At lower luminosities Lfuv decreases 
as Lir decreases, causing the simulated galaxies to lie above the 
Lfuv = 4 x 10 10 L relation. In particular, the quiescently star- 
forming galaxies tend to lie above the relation because they are 
much less obscured than the starbursts. 



3.4 SFR-M* relation 

As discussed in Section [TT] if the tight SFR-M* relation observed 
is set by the gas supply rate then galaxies that lie above the relation 
must be undergoing transient events that temporarily boost their 
SFR above what can be sustained over long time periods. Major 
mergers are one type of event that can cause galaxies to move above 
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the relation. We plot the SFR-M* relation for our simulated galax- 
ies in Fig. |11 1 along with the observed relation for z = 2.0 — 2.5 
from Karim et al. (2011 solid line). Almost all objects above the 
observed relation are starbursts, and a large fraction of the star- 
bursts are above the relation. For a given M* the starbursts can have 
SFR values > 10 x those of the quiescently star-forming galaxies. 

Note that some of the simulated galaxies, including starbursts, 
lie significantly below the relation. This is partially caused by the 
idealised setup of the simulations: Since there is no cosmological 
gas accretion some of the simulated galaxies have gas fractions 
lower than what is expected from observations and cosmological 
simulations. For fixed galaxy size and total mass, the observed KS 
relation implies that SFR scales super-linearly with gas fraction, so 
simulated galaxies with gas fractions lower than real galaxies will 
have significantly lower SFR than observed. Note also that the pre- 
burst discs have systematically higher gas fraction because the gas 
fraction decreases monotonically with time in the simulations. This 
biases the quiescently star-forming galaxies' SFR high relative to 
the starbursts, so the magnitude of the difference between starbursts 
and quiescently star-forming discs' SFR shown in Fig.[TT]should be 
taken as a lower limit. Furthermore, discrepancies in the SFR-M* 
relation of the simulations and that observed may occur because 
the SFR derived using a given diagnostic is not equivalent to the 
instantaneous SFR of the simulated galaxies, as even the observed 
relations can differ depending on what SFR diagnostic is used to 
derive them. 

Observationally, whether SMGs lie on the SFR-M* relation 
depends on the measured Lir used to infer the SFR and the in- 
ferred M*. The latter is especially difficult to determine: different 
authors have inferred masses differing by a factor of ~ 6 x for the 
same SMGs (Michalowski, Hjorth, & Watson 2010a, Michalowski, 



Watson, & Hjorth|2010b| Michalowski etal.p011| 



Hainline et al. 
201 1\ . If the Michalowski et al. (2011 ) masses are used, SMGs lie 



much closer to the SFR-M* relation than they do if the Hainline 
|et al.| ( |2011| > masses are used. However, in both cases a subset of 
SMGs are significant outliers from the relation; this conclusion is 
consistent with our claim that the SMG population is a mix of both 
quiescently star-forming and starburst galaxies. 



be problematic. For example, as discussed in Section [3^2] locally it 
seems that the CO-H 2 conversion factor Xco differs for ULIRGs 
(i.e., merger-induced starbursts) and quiescently star-forming disc 
galaxies (e.g., |Solomon et~aL]|1997[ |Downes & Solomon|[T998| ). 
If one wishes to, e.g., study possible evolution of the KS law 
with redshift then, lacking other options, it is necessary to assume 
some CO-H2 conversion factor for the high-redshift galaxy popula- 
tions observed. Choosing an appropriate Xco requires determining 
whether the high-redshift galaxies are analogous to local merger- 
induced starbursts or quiescently star-forming discs. For example, 
since it is commonly thought that SMGs are merger-driven star- 
bursts, Daddi et al. (2010) and Genzel et al. (2010) use the starburst 
Xco value for SMGs. If, however, SMGs are a mix of quiescently 
star-forming, early-stage mergers and late-stage, merger-induced 
starbursts, as we have argued in HI 1 and above, then a single Xco 
value is not appropriate for the population. In this case, use of the 
ULIRG Xco value will artificially accentuate the apparent differ- 
ences between SMGs and more typical galaxies at high redshift. 

These diagnostics can also be used to distinguish the quies- 
cently star-forming sub-populations of SMGs (blended galaxy pairs 
and isolated discs) from the starbursts in order to test the claims of 
our model and to understand the true nature of the population. We 
argue that galaxy pairs must contribute significantly to the SMG 
population (Haywar d et al.|2011b| Hay ward et al., in preparation) 
because of the weak scaling of submm flux with SFR in starbursts 
and the significantly longer duration of the galaxy-pair phase. How- 
ever, given the modelling uncertainties it is crucial to observation- 
ally determine the relative contributions. 

Furthermore, if one wishes to understand which mode of star 
formation dominates the SFR density of the universe one must be 
able to separate the modes. Even when one can clearly identify 
mergers (e.g., by the presence of tidal features) one cannot assume 
that those galaxies are dominated by merger-induced star forma- 
tion, as the SFR elevation caused by the mutual tidal torques is sig- 
nificant only near coalescence and, depending on the gas content 
and bulge fraction of the progenitors, possibly first passage. (See 
|Hopkins et al.|2010| and |Hopkins & Hernquist|20T0l for further dis- 
cussion of the distinction between star formation in mergers and 
merger-induced star formation.) The problem is amplified at higher 
redshifts when mergers cannot be easily identified. 



4 DISCUSSION 

4.1 The need to distinguish star formation modes 

If one wishes to understand star formation it is crucial to look 
beyond the local universe, because the SFR density of the uni- 
verse was greatest at z ~ 2 — 3 (e.g., |Madau et ah]|1996[|Stei 



del et al.|[T996| [Hopkins l|2004MHopkins & Beacom||2006[ |Hop^ 
kins et al.||201o| |Karim et d.||2011| |Magnelli et aL||201l| . Fur- 
thermore, the bulk of the star formation at those redshifts was ob- 
scured (e.g., Bouwens et al. 2011 , Magnelli et al. 2011), so study- 
ing IR-luminous galaxies at those redshifts is crucial. Unfortu- 
nately, galaxies become fainter and physical resolution poorer as 
one moves from z ~ to higher redshift, so observations of high- 
redshift galaxies are significantly less detailed than for local galax- 
ies. It is thus tempting to use wisdom gleaned from detailed obser- 
vations of local galaxies to guide the interpretation of observations 
of high-redshift galaxies. This is perfectly acceptable if the only 
difference between local galaxies and those at z ~ 2 — 3 is that the 
latter are farther away. However, this is clearly not the case, so one 
must apply local-universe-derived wisdom with caution. 

Assuming what is true locally is also true at z ~ 2 — 3 can 



4.2 An observational roadmap to determine what star 
formation mode powers high-redshift ULIRGs 

Fortunately, the integrated SED of a galaxy contains much infor- 
mation about the star formation mode powering it, so it is possi- 
ble to use the diagnostics we have presented here to observation- 
ally disentangle what star formation mode dominates high-redshift 
ULIRGs. This can be achieved by applying the diagnostics to the 
results of FIR and (sub)mm wide-field surveys. Since the diagnos- 
tics rely on integrated data alone they are robust to blending in 
the FIR and (sub)mm, but the beam sizes at different wavelengths 
should be similar because otherwise blending will be more severe at 
the wavelengths where resolution is poorer. The FIR and (sub)mm 
data are enough to use the Td — Lir relation as a diagnostic. Some 
of the other diagnostics require data at shorter wavelengths (to de- 
termine Lfuv and M*), so care must be taken to include all sources 
that contribute to the FIR emission, not just one component of a 
multiple component system. Keeping in mind the caveats discussed 
above, one can also estimate gas masses and use the SFE as a diag- 
nostic. 

Once ALMA SMG surveys are available it will be simple to 
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divide SMGs into single-component and multiple-component sub- 
classes. This information can be combined with Lir and T& val- 
ues (preferably derived using the PL T-distribution model) to effec- 
tively divide the SMG population into the various sub-populations: 
galaxy-pair SMGs will be resolved into two components, and those 
with smaller separations between components should typically 
have higher Lir and Td values. The merger-induced starbursts and 
isolated discs will only have one component, but the former can be 
distinguished by their higher Td values. Observations at rest-frame 
UV-NIR wavelengths can be used to apply the IRX and SFR-M* 
diagnostics to further check the classification done in the above 
manner, as starbursts will have higher values of IRX and lie above 
the SFR-M* relation. If gas masses are available one can also use 
the SFE as a diagnostic. 

Spatial information beyond the number of components can 
further aid the classification. Though IFU spectrograph data pro- 
vide high-resolution kinematics, high-redshift ULIRGs can be op- 
tically thick well into the IR, so even rest-frame near-IR observa- 
tions may not probe the central regions and thus must be interpreted 
with cautionpj Molecular gas emission, on the other hand, suffers 
significantly less dust attenuation, so (sub)mm interferometry with, 
e.g., ALMA can provide a direct view of the central regions (but 
see|Papadopoulos et al. 2010a. b). The close pairs and isolated discs 
are most effectively distinguished via interferometry because the 
former will be resolved into two components with disc-like kine- 
matics (see, e.g., Engel et al. 2010) whereas the latter will only 
show one disc. The merger-induced starbursts may show more dis- 
ordered kinematics, but it is not always simple to distinguish merg- 
ers shortly after coalescence from discs (see |Robertson & Bullock 
2008). Luckily, the starbursts can be distinguished from discs using 
the diagnostics presented here. 

By combining various data sets in the manner described 
above, it should be possible to separate high-redshift ULIRGs into 
starburst and quiescently star-forming sub-populations. In addition, 
mergers can be roughly divided into widely separated pairs and 
close pairs, which are both quiescently star-forming, and starbursts 
near coalescence. One can then more efficiently target sources for 
detailed follow-up, focusing on, e.g., different stages in the merger 
process. Furthermore, by applying this technique to SMGs one can 
compare the sizes of the sub-populations to test our claim that the 
SMG population is heterogeneous. 



4.3 Physical differences between SMGs and hot-dust 
ULIRGs 

Throughout this work we have compared the properties of qui- 
escently star-forming and starburst SMGs with simulated galax- 
ies that would not be selected by the SMG selection ('hot-dust 
ULIR Gs' ;|Chapman et al.|2004[|Chapman et al.|2008| |Casey et al.| 
|2009l|Casey et al.|2011[ |Magnelli e tal.|2010[ M12). At a given L m 
the hot-dust ULIRGs tend to have higher effective dust temperature 
than the SMGs; this is simply a consequence of the SMG selection 
because, for fixed /3 and Lir, submm flux can be decreased only by 
increasing Td. Furthermore, at a given Lir the hot-dust ULIRGs 
have higher SFE values than the SMGs, as suggested by Chapman 
|et al.| ( |2008l ), but the IRX values are similar. 



The relative locations of the SMGs and hot-dust ULIRGs on 
the SFR-M* diagram (Fig. [TTJ provide insight into the physical 
differences between these two galaxy classes: The most massive 
galaxies are almost all selected as SMGs, regardless of whether 
they are quiescently star-forming or starbursts, because they are 
luminous enough to have Ssso > 5 mJy for any reasonable Td. 
At intermediate masses almost all of the starbursts are selected as 
SMGs, but the quiescently star-forming galaxies are not. At the 
lowest masses simulated almost no galaxies are selected as SMGs. 
This reflects that fact that galaxy mass is an important driver of the 
observed submm flux because it affects both the SFR and the dust 
mass (see also Hll; Michalowski et al. 2011). The smaller galax- 
ies can be very luminous in the IR if they are undergoing a strong 
starburst, but this results in a relatively hot SED (because Td tends 
to increase sharply during starbursts; Hll) and causes them to be 
hot-dust ULIRGs rather than SMGs. 

Furthermore, Hll noted that the observed-frame submm flux, 
and thus whether a simulated starburst is selected as an SMG or 
considered a hot-dust ULIRG, depends on the angle from which 
the galaxy is viewed. Variation in the SED with viewing angle can 
cause Td to vary for different lines -of- sight, but the variation in Td 
(at most a few degrees) cannot account for the entire variation in 
submm flux. Instead, the primary cause is that the effective emis- 
sion area depends on viewing angle. Note that in Fig. [TT] there are 
many cases where large blue squares and small blue squares are 
centred at the same (M*,SFR) point. These correspond to simulated 
starbursts that are classified as both SMGs and hot-dust ULIRGs 
depending on the angle from which they are viewed. Such cases 
are relatively common, but there are no quiescently star-forming 
galaxies for which this occurs; this is most likely because in the 
quiescently star-forming galaxies most of the FIR emission arises 
from regions that are optically thin in the FIR, and thus the effec- 
tive emitting area in the FIR does not depend strongly on viewing 
angle. This interpretation is supported by the results presented in 
Section [331 



10 Rothberg & Fischer (2010 1 have demonstrated that at A = 0.65 fim dust 
obscures the nuclear discs of young stars in local (U)LIRGs; since the galax- 
ies studied here are more obscured this effect should be more severe for 
high-redshift ULIRGs. 



4.4 Other potential diagnostics 

This work has focused on diagnostics that rely on integrated UV- 
mm SEDs alone, but there are various other diagnostics that could 
potentially be used to distinguish quiescently star-forming and star- 
burst galaxies. For example, the late- stage merger-induced starburst 
SMGs may have higher CO line-widths than the quiescently star- 
forming disks (Narayanan et al. 2009). In addition, the CO spectral 
line energy distributions (SLEDs) may differ, with the starbursts 
typically having higher mean excitation than the quiescently star- 
forming galaxies. Note, however, that the significant dust optical 
depths present in the starbursts can result in SLEDs that appear in- 
dicative of very low excitation conditions, potentially making it dif- 
ficult to distinguish highly-obscured high-excitation gas from un- 
obscured low-excitation gas (Papadopou los et al.|2010a|b] >. X-ray 
observations may also be useful: In our simulations the merger- 
induced starbursts typically coincide with strong AGN activity, 
though the peak AGN activity occurs ~ 50—100 Myr after the peak 
SFR (Hopkins 2012). Thus the merger-induced starbursts should 
have stronger central X-ray sources than the quiescent discs. Fur- 
thermore, since the starbursts are typically more highly obscured, 
the AGN they host should have higher hardness ratios than the 
AGN hosted by quiescently star-forming galaxies. Finally, radio 
observations can provide excellent spatial resolution and thus be 
used to determine the number of components and separation of the 
nuclei. 



© 2012 RAS, MNRAS 000,[l]J2T] 



1 8 Hayward et al. 



5 CONCLUSIONS 

We have combined high-resolution 3-D hydrodynamic simulations 
of z ~ 2 major mergers of disc galaxies and 3-D Monte Carlo 
dust RT calculations to investigate the differences between quies- 
cently star-forming and starburst galaxies. We have focused on the 
SMG population as a case study because, as argued in Hll and 
elaborated in Section |1.2| the SMG population is likely a mix of 
quiescently star-forming galaxies and merger-induced starbursts. 
Our models make robust observational predictions for how qui- 
escently star-forming galaxies on the so-called 'main sequence of 
star formation' should differ from merger-induced starbursts. We 
present multiple observational diagnostics which can distinguish 
quiescently star-forming and starburst SMGs based on integrated 
UV-mm data alone. The testable predictions and observational di- 
agnostics presented in this work include: 

(i) Effective dust temperature - derived from fitting a single- 
T OT modified blackbody, the full form of the single-T modified 
blackbody, or a model assuming a power-law distribution of dust 
temperatures - correlates with Lir, and the galaxies in the high- 
Lir, high-Td region of the Ta — Lir plot are almost exclusively 
merger-induced starbursts. A T& cut is an effective means to select 
starbursts from the overall galaxy population, whereas an Lir cut 
only works well for those galaxies that are also selected as SMGs. 

(ii) Star formation efficiency, SFE = LiR/M gas , is ~ 2 — 3x 
higher for starburst SMGs than quiescently star-forming SMGs at 
a given Lir. The starbursts not detected as SMGs have SFE ~ 2x 
greater than the starburst SMGs. 

(iii) The IR excess, IRX = Lir/Lfuv, correlates with Lir. 
The starbursts have IRX a factor ~ 5x greater than quiescently 
star-forming galaxies, primarily because they are typically more 
luminous than quiescently star-forming galaxies and IRX correlates 
with Lir. At the bright end the correlation is approximately linear, 
suggesting that Lfuv is decoupled from Lir for the most obscured 
starbursts. 

(iv) Effective dust temperature, star formation efficiency, and 
IRX are all inversely correlated with nuclear separation because 
the strength of the merger-driven starburst increases as the galaxies 
coalesce. 

(v) The majority of the simulated starbursts lie above the SFR- 
M* relation, whereas the quiescently star-forming galaxies lie close 
to it. The observational definition of starbursts as outliers above 
the relation agrees well with our definition of merger-induced star- 
bursts as simulated galaxies with instantaneous SFR > 3x their 
baseline quiescent SFR. 

One can apply these observational diagnostics to test our claim 
that the SMG population is a mix of quiescently star-forming galax- 
ies and merger-induced starbursts and to constrain the relative con- 
tribution of the sub-populations. Furthermore, the tests presented 
here provide physically motivated ways to observationally separate 
quiescently star-forming galaxies and starbursts, enabling one to 
more cleanly study the underlying physics than when heterogenous 
samples (e.g., SMGs) are used. Though we have focused on the 
SMG population here, we have also discussed the hot-dust ULIRGs 
that would not be selected by the SMG selection. These galaxies 
tend to be less massive than SMGs. Though they can be as lumi- 
nous as SMGs, this often requires that they be starbursts, so they 
have relatively high Td for their Lir and are thus missed by the 
SMG selection. Hot-dust ULIRGs are also sometimes just SMGs 
viewed from a different angle, because the submm flux of a simu- 
lated starburst can vary significantly with viewing angle. 



We have also explored how well various IR SED fitting forms 
describe our simulated galaxies' SEDs. For a given SED, the fitted 
Td and f3 can vary significantly depending on what form is used, 
so these parameters should not be interpreted physically. How- 
ever, some trends are robust to choice of fitting method, indicating 
that the fits still yield useful information about the SED variations 
amongst galaxies. In future work we will more thoroughly investi- 
gate how well the different FIR SED fitting methods describe ob- 
served and simulated galaxy SEDs, determine how the inferred dust 
properties relate to the actual physical dust properties, and develop 
an improved method for fitting the FIR SEDs of large samples of 
galaxies. 
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